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Di-jet production is studied in collisions of quasi-real photons at e"'"e~ centre-of-mass energies y/Sce from 189 
to 209 GeV at LEP. The data were collected with the OPAL detector. The structure of jets is investigated and 
differential cross sections are measured and compared to QCD calculations. 



1. Introduction 

We have studied the production of di-jets in the 
collisions of two quasi-real photons at an e^'e^ 
centre-of-mass energy y^s^ from 189 to 209 GeV, 
with a total integrated luminosity of 593 pb~^ 
collected by the OPAL detector at LEP. Di-jet 
events are of particular interest, as the two jets 
can be used to estimate the fraction of the photon 
momentum participating in the hard interaction, 
which is a sensitive probe of the structure of the 
photon. The transverse energy of the jets pro- 
vides a hard scale that allows such processes to be 
calculated in perturbative QCD. Fixed order cal- 
culations at next-to-leading order (NLO) in the 
strong coupling constant as for di-jet production 
are available and are compared to the data, pro- 
viding tests of the theory. Leading order Monte 
Carlo (MC) generators are used to estimate the 
importance of soft processes not included in the 
NLO calculation. More details on the results pre- 
sented here can be found in P . 

The fc_L-clustering algorithm [2] is used as op- 
posed to the cone algorithm |3] in our previous 
publications |4I5| for the measurement of the dif- 
ferential cross-sections, because of the advantages 
of this algorithm in comparing to theoretical cal- 
culations inj. The cone jet algorithm is used to 
demonstrate the different structure of the cone 
jets compared to jets defined by the fcj^ -clustering 
algorithm. 

At e+e" colliders the photons are emitted 
by the beam electrons (positrons). Most of 
these photons carry only a small negative four- 
momentum squared, Q^, and can be considered 



quasi-real (Q^ « 0). The electrons are hence scat- 
tered with very small polar angles and are not 
detected. Events where one or both scattered 
electrons are detected are not considered in the 
present analysis. 

In LO QCD, neglecting multiple parton inter- 
actions, two hard parton jets are produced in 77 
interactions. In single- or double-resolved interac- 
tions, these jets are expected to be accompanied 
by one or two remnant jets. A pair of variables, 
xt and be defined J that estimate the 

fraction of the photon's momentum participating 
in the hard scattering: 



_|_ jets— 1,2 

hfs 



(1) 



where Pz is the momentum component along the 
z axis of the detector and E is the energy of 
the jets or objects of the hadronic final state 
(hfs). In LO, for direct events, all energy of the 
event is contained in two jets, i.e., — 1 and 
x~ — I, whereas for single-resolved or double- 
resolved events one or both values are smaller 
than 1. Differential cross sections as a function 
of x.y or in regions of Xj are therefore a sensitive 
probe of the structure of the photon. 

2. Jet structure 

The internal structure of jets is studied using 
the jet shape, which is defined as the fractional 
transverse jet energy contained in a subcone of 
radius r concentric with the jet axis, averaged 
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Figure 1. The jet shape, '^{r), for the two regions 
of a;+-x~-space indicated in the figure (a), and 
^'(r) for < 0.75 compared to the predictions 
of the LO MC generators PHOJET and PYTHIA 
(b). Figures (c) and (d) show the value of 'i'{r — 
0.4) as a function of the transverse energy and 
pseudo-rapidity of the jet respectively, compared 
to the PYTHIA prediction. 



Figure 2. The value of the jet shape \E'(r) at 
r = 0.4 as a function of the jet transverse energy 
for < 0.75 (a) and > 0.75 (b), and as a 
function of the jet pseudo-rapidity for x^ < 0.75 
(c) and x^ > 0.75 (d). In each figure the re- 
sults obtained using the inclusive k± and the cone 
jet algorithm are shown and compared to the 
PYTHIA prediction. 



over all jets of the event sample: 

^ B¥'(r = LO) 

with r = a/(A7?)2 -I- (A0)2 and A'jot the total 
number of jets analysed. Both k±_ and cone jets 
are analysed in this way. As proposed in [HI , only 
particles assigned to the jet by the jet finders are 
considered. Events entering the jet shape distri- 
butions are required to have at least two jets with 
a transverse energy 3 GeV < Ei^^ < 20 GeV and 
a pseudo-rapidity \ri^°^\ < 2. 

In Figure ni^) the jet shape, ^'(r), is shown 
for the kj_ algorithm for both x^ > 0.75 and 
x^ < 0.75. Here and in subsequent figures the 
total of statistical and systematic uncertainties 
added in quadrature is shown where larger than 
the marker size. The inner error bars show the 



statistical errors. The first sample is dominated 
by direct photon-photon interactions and hence 
by quark-initiated jets. As is demonstrated in 
the figure, jets in this sample are more coUi- 
mated than for small values of x^ , where the 
cross-section is dominated by resolved processes 
and hence has a large contribution from gluon- 
initiated jets. In both cases the jets become more 
collimated with increasing transverse energy, as 
is shown in Figure ^(c). There is no signifi- 
cant dependence on the jet pseudo-rapidity (Fig- 
ure [Tl(d)). Both PHOJET and PYTHIA ^ 
give an adequate description of the jet shapes as 
can be seen in Figures n](b) , (c), and (d). The 
default choices of SaS ID [TT] for PYTHIA and 
LO GRV 12 for PHOJET are taken. 

Figure |21 compares the shapes of jets defined 
by the cone algorithm and the k±^ algorithm, in 
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each case compared to the shape as obtained from 
PYTHIA. As for the A: -jets, the jets defined 
by the cone algorithm are more coUimated in 
the quark-dominated sample and always become 
more collimated for increasing transverse energy, 
while there is no dependence on the jet pseudo- 
rapidity. The cone-jets are significantly broader 
than the jets defined by the kj_ algorithm at low 
Ei^\ With increasing Ei^*, jets become more col- 
limated and the two jet algorithms give similar 
results. While the fcj^-jets are well described by 
PYTHIA and PHOJET, the jet shapes obtained 
for the cone-jets are somewhat broader than in 
the data. 

3. Differential Di-jet cross-sections 

Only the k±_ jet algorithm is used for the mea- 
surement of the differential di-jet cross-sections. 
The experimental results are compared to a per- 
turbative QCD calculation at NLO ^31 which 
uses the GRVHO parametrisation of the par- 
ton distribution functions of the photon ^21 : and 
was repeated for the kinematic conditions of the 
present analysis. The renormalisation and fac- 
torisation scales are set to the maximum Ei^^ 
in the event. The calculation was performed 
in the MS-scheme with five light flavours and 
A^^D = 130 MeV. The average of the hadronisa- 
tion corrections estimated by PYTHIA and HER- 
WIG have been applied to the calculation for this 
comparison. In the figures described below the 
shaded band indicates the theoretical uncertainty 
estimated by the quadratic sum of two contribu- 
tions: variation of the renormalisation scale by 
factors of 0.5 and 2 and the difference between 
using HERWIG or PYTHIA in estimating the 
hadronisation corrections. 

Due to the different nature of the underlying 
partonic process one expects different distribu- 
tions of the angle Q* between the jet axis and 
the axis of the incoming partons or direct pho- 
tons in the di-jet centre-of-mass frame. The lead- 
ing order direct process 77 qq proceeds via the 
t-channel exchange of a spin-i quark, which leads 

to an angular dependence oc (l — cos^ e*) . In 
double resolved processes the sum of all matrix el- 
ements, including a large contribution from spin-I 
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Figure 3. The di-jet cross-section as a function of 
|cos9*| for the two regions in a;+-a;~-space indi- 
cated in the figure. 



gluon exchange, leads to an approximate angular 
dependence oc (1 — |cos6*|) ^ J^l- The contri- 
bution of the different processes to all resolved 
events depends on the parton distribution func- 
tions of the photon. An estimator of the angle O* 
can be formed from the pseudo-rapidities of the 
two jets as 

(ict jet \ 

'"^^ , (3) 

where it is assumed that the jets are collinear in (/) 
and have equal transverse energy. Only |cosS*| 
can be measured, as the ordering of the jets in 
the detector is arbitrary. To obtain an unbiased 
distribution of |cosO*| the measurement needs to 
be restricted to the region where the di-jet in- 
variant mass Mjj = 2El^^ j 1 — |cos 9*|^ is not 

influenced by the cuts on E}^^ In the present 
analysis a cut of Mjj > 15 GeV ensures that 
the |cos6*| distribution is not biased by the re- 
strictions on E}^ for the range |cosO*|<0.8 and 
|7^j°t| = I ^77!°' + ?72°') /2| < 1 confines the mea- 
surement to the region where the detector resolu- 
tion on |cos8*| is good. 

Figure |21 shows the differential di-jet cross- 
section as a function of |cos6*| for both > 
0.75 and < 0.75. The steeper rise with 
increasing |cosO*| from the dominating spin-1 
gluon exchange in the second sample is clearly vis- 
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Figure 4. The di-jet cross-section as a function 
of the mean transverse energy E^r^^ of the di-jet 
system, for the three regions in a;;:Jr-a;~-space given 
in the figure. The factor / is used to separate the 
three measurements in the figure more clearly. 



Figure 5. The di-jet cross-section as a function 
of x-y and for the regions of the mean transverse 
energy Ei^^ and of the di-jet system indicated 
in the figures. 



ible. The shape of both samples is well described 
by NLO QCD. For x^ < 0.75 the NLO calcula- 
tion is about 20% below the data. It should be 
noted that in this region the contribution from 
the underlying event, not included in the cal- 
culation, is expected to be largest, as discussed 
m more detail below. For x^ > 0.75 the NLO 
QCD prediction is about 20% above the data. 
While here the contribution from MIA is small, 
this region is affected by rather large hadronisa- 
tion corrections, which translates into an uncer- 
tainty of the normalisation in comparing the the- 
oretical prediction to the data. Probably more 
importantly it has been pointed out that the cal- 
culation of the cross section becomes increasingly 
problematic when approaching x^ = 1 jl5j . 

The differential di-jet cross-section as a func- 
tion of the mean transverse energy Ei^^ of the 
di-jet system is shown in Figure H At high e)^^ 
the cross-section is expected to be dominated 
by direct processes, associated with the region 



x^ > 0.75. Consequently we observe a signif- 
icantly softer spectrum for the case x^ < 0.75 
than for the full x+-a;~-space. The calculation 
is in good agreement with the data for the full 
x^ -x~ -range and for x^ or x~ < 0.75. The cross- 
section predicted for x^ < 0.75 is again below the 
measurement. PYTHIA 6.161 is in good agree- 
ment with the measured distributions using the 
SaS ID parton densities. 

The three plots of Figure El show the differ- 
ential cross section as a function of x.y for the 
three regions in xj^-x~-spa.ce described above. 
The shaded histogram on the bottom of each 
of the three plots indicates the contribution of 
MIA to the cross section as obtained from the 
PYTHIA MC generator. It is evident espe- 
cially for x^<l that the MIA contribution is of 
about the same size as the discrepancy between 
the measurement and the NLO prediction. Fur- 
thermore it is interesting to observe that there is 
next to no MIA contribution to the cross section 
if either x^ or x~ is required to be less than one. 
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while the sensitivity to the photon structure at 
small is retained. As one would expect also 
the agreement of the NLO calculation with the 
measurement is best in this case. With these 
measurements one is therefore able to disentangle 
the hard subprocess from soft contributions and 
make the firm statement that NLO perturbative 
QCD is adequate to describe di-jet production in 
photon-photon collisions in the regions of phase 
space where the calculation can be expected to be 
complete and reliable, i.e. where MIA contribu- 
tions are small and for not too close to unity. 
At the same time a different sub-set of observables 
can be used to study in more detail the nature of 
the soft processes leading to the underlying event. 
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